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We introduce a family of 2D topological subsystem quantum error-correcting codes. The gauge 
group is generated by 2-local Pauli operators, so that 2-local measurements are enough to recover 
the error syndrome. We study the computational power of code deformation in these codes, and 
show that boundaries cannot be introduced in the usual way. In addition, we give a general mapping 
connecting suitable classical statistical mechanical models to optimal error correction in subsystem 
stabilizer codes that suffer from depolarizing noise. 
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I. INTRODUCTION 

Quantum error correction [l|-[3l and fault-tolerant 
quantum computation promise to allow almost per- 
fect storage, transmission and manipulation of quantum 
information. Without them, quantum information pro- 
cessing would be doomed to failure due to the decoher- 
ence produced by interactions with the environment and 
the unavoidable inaccuracies of quantum operations. 

The key concept in quantum error correction is the no- 
tion of quantum code. This is a subspace of a given quan- 
tum system where quantum information can be safely en- 
coded, in the sense that the adverse effects of noise can be 
erased through an error correction procedure. In practice 
this procedure is also subject to errors and thus it should 
be as simple as possible to minimize them. Naturally, 
the meaning of 'simple' will depend on particular imple- 
mentations. A common situation is that interactions are 
restricted to quantum subsystems that are close to each 
other in space. In those cases, the locality of the opera- 
tions involved in error correction becomes crucial. 

The stabilizer formalism [Tol . ITlT ] provides a unified 
framework for many quantum codes. In stabilizer codes 
the main step for error correction is the measurement of 
certain operators, which may be local or not. A class of 
codes where these measurements are intrinsically local is 
that of topological stabilizer codes [l2l - [l5| . In a different 
direction, locality can also be enhanced b y c onsidering 
more generally stabilizer subsystem codes [III li~7j . The 
present work provides an example of a family of codes 
which can be labeled both as 'topological' and 'subsys- 
tem'. 

Topological codes where originally introduced with the 
goa l of obtaining a self-protecting quantum computer 
12]. This idea faces important difficulties in low di- 
mensions, since thermal instabilities are known to occur 
EE El- On the other hand, topological codes are 
local in a natural way and have very interesting features 
in the context of active error correction. For example, 
they do not only allow to perform operations transver- 
sally [l3, HU but also through code deformations [H, [2(| ■ 
Moreover, there exist a useful connection between error 
correction in topological codes and certain classical sta- 



tistical models [la, |21j, |22j . 

Stabilizer subsystem codes are the result of applying 
the stabilizer formalism to operator quantum error cor- 
rection [la] . In subsystem codes part of the logical qubits 
that form the code subspace are no longer considered as 
such but, rather, as gauge qubits where no information 
is encoded. This not only allows the gauge qubits to ab- 
sorb the effect of errors, but has interesting consequences 
for error correction. It may allow to break up each of 
the needed measurements in several ones that involve a 
smaller number of qubits Uv i [l 7ll . An example of this is 
offered by Bacon-Shor codes 16] , in which the basic oper- 
ators to be measured can have support on an arbitrarily 
large number of qubits, yet their eigenvalues can be re- 
covered from 2-local measurements that do not damage 
encoded information. Moreover, the pairs of qubits to 
be measured together are always neighbors in a 2D lat- 
tice. Thus, subsystem codes can have very nice locality 
properties. 

The 2D topological subsystem codes introduced here 
show all the characteristic properties of topological codes 
and at the same time take profit of the advantages of 
subsystem codes. Some of them are: 

• The codes are local in a 2D setting, which can be 
flat. 

• The measurements needed for error correction only 
involve two neighboring qubits at a time, as in 
Bacon-Shor codes. This is an important advan- 
tage with respect to other topological codes, such 
as surface codes, that require measuring groups of 
at least 4 qubits. 

• Most errors of length up to an are correctable, 
where c is some constant and n the number of phys- 
ical qubits. This feature, common to topological 
codes, follows from the fact that logical operators 
are related to strings that wind nontrivially around 
the surface where the code is defined. 

• Error correction must be done only 'up to homol- 
ogy', an important simplification that allows the 
introduction of specific tools. 



2 



• One can naturally perform certain logical gates 
through 'deformations' of the code. This feature, 
however, is less powerful that in other topological 
codes because boundaries cannot be introduced in 
the usual way. 

Since these codes are topological, it is natural to ex- 
pect a connection between their correction procedures 
and suitable classical statistical models. However, the 
ma ppi ng between surface codes and random Ising mod- 
els [13} makes stronge use of their CSS structure [3, [25[ , 
and the same is true for the one between color codes 
and 3-body random Ising models |22j. The CSS struc- 
ture makes it possible to completely separate the correc- 
tion of phase flip and bit flip errors, making the problem 
classical and enabling the connection. Indeed, there ex- 
ist similar mappings from classical codes to statistical 
models (26|. Fortunately, as explained in section VI, the 
approach can be generalized even in the absence of this 
separation. Moreover, the subsystem structure is also 
compatible with the approach, so that it can be applied 
to the family of codes of interest. 

The paper is organized as follows. Sections II and III 
go over several aspects of quantum error correction and 
topological codes, respectively, setting up a framework 
for the rest of the paper. Section IV introduces the family 
of topological subsystem codes and presents a thorough 
study of their properties. Section V offers the construc- 
tion of a general mapping between error correction in 
subsystem codes and classical statistical models. Section 
VI is devoted to conclusions. 



II. STABILIZER QUANTUM 
ERROR-CORRECTING CODES 

This section summarizes the notions of quantum error 
correction that will be needed in the rest of the paper. 
It mainly reviews stabilizer codes, both in the subspace 
and the more general subsystem formulation. Ideal error 
correction procedures and their success probability are 
also considered. 



A. Quantum error correction 

Quantum error correction deals with the preservation 
of quantum information in the presence of noise. Both 
the noise £ and the error recovery TZ are modeled as 
quantum operations or channels £,TZ : L(H) — > L(H), 
where L(H) is the space of linear operators on H, the 
Hilbert space associated to the quantum system under 
consideration. Such maps can always be expressed in the 
operator-sum representation. For example, the noise is 
£{p) = J^i^iP^i f° r somc Ei e L(H), which will be 
denoted by £ — {Ei}. 

In the original formulation of quantum error correction 
(SHU j quantum information is encoded in a subspace of 



H, the code subspace C c H. The system undergoes 
a noisy process £ and afterwards an error recovery op- 
eration TZ is performed. Then, given a code C, a noise 
source £ is said to be correctable if there exists a recov- 
ery operation TZ such that 1Z o £{p) = p for any state 
P eL(C). 

More generally, in the operator quantum error correc- 
tion formalism [231 ]. information is encoded in a subsys- 
tem A, with C = A&B. Whatever happens to subsystem 
B is irrelevant. That is, error recovery is possible for a 
quantum channel £ if there exists a recovery operation 
TZ such that for any p A G L(A) and p B e L(B) it gives 
TZ o £(p A ® p B ) — p A ® p' B for some arbitrary p' B . 

The sufficient and necessary condition for the noise 
process £ = {Ei} to be correctable [1, 0, [23| is that 
PEjEjP = l A ®g B for every i and j, with P the projec- 
tor onto the code subspace. When this condition holds, 
the set of errors {Ei} is said to be correctable. Since 
adding a linear combination of the E a to the set does 
not change correctability, it is natural to consider cor- 
rectable sets of errors as linear subspaces and to choose 
the most convenient operator basis. Generally the quan- 
tum system is composed of n qubits, H ~ (C 2 )® n , and 
error operators are chosen to be Pauli operators, elements 
of the Pauli group V n := (il, X\, Z\, . . . , X n , Z n ). Here 
Xi, Zi are as usual the Pauli operators on the i-th qubit, 
X= |0)(1| + |1}(0|, Z= |0> <0| — 1 1> <1 1 in the orthonormal 
basis {|0), |1)}. 

Usually error models are such that errors which affect 
more qubits are less likely to happen. Then it makes 
sense to correct as many errors as possible among those 
that have support on (act nontrivially on) a smaller num- 
ber of qubits. The weight \E\ of a Pauli operator E € V n 
is defined as the number of qubits that form its support. 
When a code can correct all Pauli errors E with \E\ < r 
it is said to correct r errors. 



B. Stabilizer subspace codes 

A formalism that has been particularly successful for 
the development of quantum codes is the stabilizer for- 
malism [13, [Tl|, m which the code C is described in terms 
of an Abelian subgroup Q = (Gj) C V n such that —1^5. 
Take the generators Sj to be independent and let s be 
the rank of S. The n qubit Hilbert space H can be parti- 
tioned according to the eigenvalues of the Sj into 2 s iso- 
morphic orthogonal subspaces H — @ s C s . Here s = (sj) 
is the error syndrome, with Sj = ±1 the eigenvalue of Sj. 
By convention, the code subspace C is that with Sj = 1 
for all j. It has dimension 2 k , with k = n — s the number 
of encoded or logical qubits. The reason to call s the er- 
ror syndrome is that in can be obtained by measuring the 
Sj and then used to infer which errors have produced. 

It is easy to introduce a Pauli group for the k logical 
qubits. Let N(S) be the normalizer of S in V n . Its 
elements are the Pauli operators that map the subspaces 
C s onto themselves, and the quotient group N(S)/S is 
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isomorphic to Vk- The logical Pauli operators are then 
generated by X 1; Z 1; . . . , X k , Z k £ N(S), some chosen 
representatives of the images of Xx, Z\, . . . , X k , Z k £ Vk 
under a given isomorphism. 

It is also possible to characterize a stabilizer code with 
a pair (U, s), where U is an automorphism of V n and s an 
integer, < s < n. Let Xi, Zi denote the images of Xj, 
Z{ trough U. Then the stabilizer is S = (Zx, • • • , Z s ). 
This approach directly provides a choice for encoded 
Pauli operators: X\ := X s+ i,Zi := Z s +x> ■ ■ ■ >Xk '■= 
An, Zk '■= Z n . 

Pauli errors E have a specially simple effect on the 
encoded states, as they map the subspaces C s one onto 
another. Set = s when ESjE^ — SjSj. Then a Pauli 
error E maps C onto C Se . Pauli operators are divided 
in three categories. The elements of V n — N(S) map the 
code to other subspaces and are termed detectable errors, 
as their effect can be detected by measuring the operators 
{Sj}. The elements of N(S) -S', with S' := (il)«S, map 
in a nontrivial way the code to itself and are thus called 
undetectable errors. Finally, the elements of S' have no 
effect on encoded states p £ L(C). The distance d of 
the code is defined as the minimum weight among un- 
detectable errors. It determines the number of corrected 
errors, which is [(d — 1)/2J . A code of n qubits that 
encodes k qubits and has distance d is denoted [[n, k, d]]. 

C. Stabilizer subsystem codes 

The stabilizer formalism can also be used in the context 
of operator quantum error correction (l6l . [l7| . Instead 
of being characterized by a stabilizer group, subsystem 
stabilizer codes are almost determined by a subgroup Q C 
V n , called the gauge group, such that il £ Q. Almost 
because, in addition, a stabilizer group S has to be chosen 
such that S', as defined above, is the center of Q. There 
are different choices for S because the sign of some of 
its generators can always be flipped. This amounts to 
different choices for C in the decomposition H = @ s C s . 

The idea after the introduction of the gauge group is 
that gauge operations should not affect the encoded in- 
formation. This forces to identify states such as p and 
GjpGj as equivalent, giving rise to a subsystem structure 
C s = A s (g) B s . The decomposition is such that the gauge 
operators Gj act trivially in the A s subsystems and gen- 
erate the full algebra of operators of the B s subsystems. 
Set C = A (g) B, with A the logical subsystem where in- 
formation is encoded and B the gauge subsystem that 
absorbs the effect of gauge operations. Since Q /S ~ V r 
for some r, B consists of r qubits. Similarly, the Pauli 
operators for the k logical qubits arc recovered from the 
isomorphism N(Q)/S ~ Vk, and k + r + s = n. 

It is also possible to characterize a stabilizer sub- 
system code as a triplet (U,s,r), where U is an au- 
tomorphism of V n and s,r > are integers with 
r + s < n. Using the same notation as above, 
the stabilizer and gauge groups are S — (Zx, ■ ■ ■ ,Z S ) 



and G = (il,Zx,...,Z s+r ,X s+ x,...,X r ). The cho- 
sen logical Pauli operators are X\ := X s+r+ i, Z\ := 
Z s +r+i, ■ ■ ■ , Xk ■— X n , Zk := Z n . 

In subsystem codes detectable Pauli errors are the el- 
ements of V n — N(S) and undetectable ones are those 
in N(S) — Q. Undetectable errors are directly related 
to logical Pauli operators. Indeed, N(S)/G ~ N(Q)/S', 
through the following correspondence. For any E £ N(S) 
there exists a G £ Q such that EG £ N(Q), and if G" £ Q 
is such that EG £ N(Q) then GG' £ S' = GnN(G). The 
distance d of the code is defined as for subspace codes and 
has the same implications regarding error correction. A 
subsystem code of n qubits that encodes k qubits and 
has r gauge qubits and distance d is denoted [[n, k, r, d]}. 

D. Syndrome measurement 

An interesting property of stabilizer subsystem codes 
is that they may allow an easier measurement of the sta- 
bilizer generators. This is so because it is possible to 
substitute the direct measurement of an stabilizer ele- 
ment S by an indirect one, in which t self-adjoint gauge 
operators Gi such that S = G\ ■ ■ ■ Gt are measured. It 
may be the case that the Gi have a smaller weight than 
S. For example, in the family of Bacon-Shor codes the 
gauge generators Gi have always weight \G\ = 2 but 
the smallest stabilizer generators have an arbitrarily large 
weight. Such cases offer two important advantages. On 
the one hand, the smaller the weight of a Pauli operator, 
the simpler the operations needed to measure it. This is 
specially relevant in fault-tolerant quantum computing, 
where error correction is considered a faulty process in 
itself, because simpler operations imply less errors. Sec- 
ondly, it may be possible to measure the Gi in parallel, 
with the corresponding saving of time. This is again rel- 
evant for fault tolerance, where the ability to perform 
measurements faster entails less errors. 

Since the Gi need not commute the ordering of the 
measurements is relevant. In general, the ordered mea- 
surement of a collection of t operators Ex , . . . , E t £ V n 
yields the effective measurement of an abelian group of 
self-adjoint Pauli operators Ai C Z, with Z the center 
of (-l,Ex,...,E t ). In particular, M = Af n Z with 
Af the abelian group of those self-adjoint Pauli opera- 
tors with eigenvalues fixed by the sequence of measure- 
ments. Af can be computed iteratively, since adding 
an additional measurement E t +\ changes Af to Af 1 = 
{E t+1 )-(AfnN(E t+1 )) 

E. Error correction 

Even if a noisy channels £ is not correctable for a given 
subsystem code, there exists some probability of perform- 
ing a successful error recovery. For example, if each of 
the physical qubits that compose a code undergoes a de- 
polarizing channel {(1 — 3p)? l,p? X,pzY,pz Z}, then the 
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noise is certainly not correctable, but the success prob- 
ability can still be close to one. This section quantifies 
this probability, which is of primary importance for topo- 
logical codes (l3l |. 

Some notation is needed here. Set as equivalent E ~ 
E' those operators E,E' E L(H) that have the same 
action up to gauge elements, E — E'G for some G € Q. 
The corresponding equivalence classes will be denoted E. 
Let {Ei\f^ x C N(S) be a particular set of representatives 
for N(S)/G, taking in particular D\ = 1. The Di with 
i > 1 will represent the ways in which error correction 
can fail. For example, if there is one encoded qubit a 
choice is {1,X\,Z\,X-lZ{\. To exted the equivalence of 
operators to channels, we choose the minimal equivalence 
relation such that {£";} ~ {E'J whenever E; k ~ E[ for all 
i. 

Assume an error model in which Pauli errors E € V n 
occur with a given probability p(E). That is, £ = 
{p(E)z E}EeVn- Errors with different phases will not 
be distinguished when discussing p(E) because phases 
are irrelevant. Up to gauge operations the error channel 
is £ ~ {p(E)z E}^ £ -p n /g, where E denotes any chosen 

element of E and p(E) = YlGeg/(ii)P(^^) * s ^ ne Prob- 
ability for a given class of errors to happen. This makes 
already apparent that class probabilities p(E) are more 
important than individual error probabilities p(E) . 

Error recovery starts with the measurement of the sta- 
bilizer generators Sj. This yields the error syndrome s, 
which limits the errors E that have possibly happened 
to those with s E = s. These possible errors are arranged 
into different classes, which may be labeled by choosing 
any possible error E and taking as representatives the el- 
ements {EDi}. Then the conditional probability for the 
class of errors E to have happened given the syndrome 
outcome s is 



p(E\s) = 



p(E) 



(i) 



Suppose that these conditional probabilities can be com- 
puted, which may be potentially difficult due to the 
combinatorics. Then the class E = E s that maximizes 
p(E\s) is known and the optimal recovery operation is 
1Z = {E S P S } S , where P s is the projector onto the sub- 
space C s . The combined effect of errors and recovery is 
i 

Tl o £ ~ {p?Di}i for some probabilities pi that only de- 
pend on the error distribution p(E). This gives a success 
probability for the error recovery 
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(2) 



A bad feature of the expression ([2]) is that it depends 
on E s . Consider an alternative error correction proce- 
dure where the class of errors with the maximum prob- 
ability E s is not always chosen. Instead, an operator 
from a class E is applied with probability p(E\s), so that 
TV = {p(E\s) z EP S } S E . The success probability for this 



randomized correction procedure is 

p' = Y,P(E)p(E\s E ) = ^2p(E) P (E\se). (3) 

E E 

This procedure is at best as successful as the original one, 
giving the bound p' < pa- Notice that po = I if and only 
if p' = 1 . It follows that po = 1 if and only if for any 
D € N(S) - g 



Y,P(E)p(ED\s E ) = 0. 



(4) 



E 



This was the condition used in [13| to characterize suc- 
cessful recovery. 



III. TOPOLOGICAL STABILIZER CODES 

This section gathers together several aspects of topo- 
logical stabilizer codes to provide a reference for section 
rvl The goal is to put the subsystem codes introduced 
there in a broader context, making apparent the simi- 
larities and differences with previously known local and 
topological codes. 



A. Local codes 

In the context of fault-tolerant quantum computing 
it is advantageous to be able to perform the syndrome 
measurements in a simple way. It may be desirable, for 
example, that the number of qubits that form the support 
of the operators to be measured is small. Similarly, it 
may be convenient that the qubits that form the code 
only belong to the support of a small number of such 
operators. This ideas can be formalized to give rise to 
the notion of local families of codes [l2]. In particular, 
a family of stabilizer subspace codes {C;} is local when 
i/ it contains codes of arbitrary distance and ii/ there 
exist two positive integers /i, v such that for each Ci there 
exists a family of generators of the stabilizer Sj such that 
ii.a/ \Sj \ < [i and ii.b/ the number of Sj-s with support 
on any given qubit is smaller than v. 

The drawback of such an abstract notion of locality is 
that it is not related in any way to a particular geome- 
try. Many physical settings have qubits disposed in 1,2 
or 3 spatial dimensions and only allow the direct inter- 
action of nearby qubits. To reflect this fact, and without 
loss of generality regarding the nature of the lattices, the 
definition above may be modified as follows [13] ■ First, 
the qubits of each code Ci are considered to be disposed 
in the vertices of some finite, periodic, cubic lattice of 
a given dimension D. Second, instead of ii/ there must 
exist a positive number d, independent of Ci, such that 
the support of all the Sj is contained in some cube con- 
taining d D vertices. A family of codes that is local in D 
dimensions is always local in the previous sense. 
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FIG. 1: In this figure the geometry is that of a torus, with 
opposite sides of the square identified. The dark region can 
be contained in a disconnected collection of simply connected 
regions, and thus cannot be the support of a non-detectable 
error in a topological code. 



A similar notion of local stabilizer subsystem codes 
may be defined by substituting the stabilizer with the 
gauge group. Notice that a family of subsystem codes 
might have local gauge generators but non-local stabilizer 
generators, as Bacon-Shor codes exemplify. 



B. Topological codes 

Topological stabilizer codes are constructed from lat- 
tices on a given manifold, in such a way that the genera- 
tors of the stabilizer are local with respect to the lattice 
[II III EE SI- Typically the set of qubits and the set of 
generators are directly related to sets of geometric objects 
such as the vertices, links or faces of the lattice. In order 
to distinguish truly topological codes from merely local 
ones, we propose the following criterium. In a topologi- 
cal stabilizer code, any operator O € N(S) with support 
in a subset of a region composed of disconnected pieces, 
each of them simply connected, has trivial action on logi- 
cal qubits. Stated this way, it is a rather vague criterium 
since no formal definition of region or connectedness is 
given. However, it will be enough for our purposes by 
adopting a reasonable interpretation when needed. Fig. [I] 
shows an example of a region in a torus that cannot be 
the support of an undetectable error. 

An enumeration of the common properties of known 
topological codes will be useful. First, the support of un- 
detectable errors is topologically non-trivial in some well- 
defined way. Indeed, both in surface codes [lU, 0, [28| 
and color codes [HI, HH undetectable errors are related 
to homologically non-trivial cycles. Second, and closely 
related, the number of encoded qubits k depends only 
upon the manifold in which the lattice is defined, and 
not the lattice itself. For example, for 2D surface and 
color codes the number of logical qubits is respectively 
k = 2 — x and k — 4 — 2%, with x the Euler characteris- 
tic of the 2-manifold. Finally, an important property of 



topological codes is that their nature makes it possible to 
define them in many different lattices, typically rather ar- 
bitrary ones as long as some code dependent constraints 
are satisfied. In other words, topological codes display 
a huge flexibility. This is to be expected from construc- 
tions that only see the topology, not the geometry, of the 
manifolds to which they are related. 

Notice that the locality of the stabilizer generators has 
been emphasized, with no mention to gauge generators. 
The reason is that, up to now, no genuinely subsystem 
topological codes have been known. This paper intro- 
duces a family of such codes. They have both local gauge 
generators and local stabilizer generators. Such locality 
properties should be expected from any topological sub- 
system code. 

The family of Bacon-Shor codes [l6| provides an ex- 
ample of non-topological local gauge codes. These codes 
certainly do not satisfy the above criterium for any in- 
terpretation of connectedness that agrees with their 2D 
lattice geometry. Moreover, their geometry is completely 
rigid, in the sense that there is no clear way to generalize 
them to other lattices and manifolds. 

It is interesting to observe that topological codes do not 
offer good d/n ratios, which go to zero as larger codes are 
considered. For example, in two-dimensional surface or 
color codes d = 0{^/n) holds (which is optimal among 
2-dimensional local codes [13]). But, as remarked in fl2^ . 
this is a misleading point because topological codes can 
correct most errors of order 0(n). 

Finally, classical topological codes also exist [28]. Un- 
like quantum ones, they can be obtained from mere 
graphs, that is, 1-dimensional objects. 



C. Topological quantum memories 

In [l3| an interesting approach to the problem of indef- 
inite preservation of quantum information was presented 
that makes use of topological codes. Since it will underly 
several discussions below, a brief summary is in order. 
The main idea is that information is encoded in a surface 
code and preserved by keeping track of errors. To this 
end, round after round of syndrome extractions must be 
performed. There are thus two sources of errors, since not 
only the code will suffer from storage errors but the stabi- 
lizer measurements themselves are also faulty. When the 
error rate is below certain threshold, the storage time can 
be made as long as desired by making the code larger, 
a feature that is only available in topological codes (for 
other codes one would have to use concatenation) . Inter- 
estingly, this error threshold can be connected to a phase 
transition in a classical statistical model, a random 3D 
gauge model [H, [H[ . 
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D. Code deformation 

The flexibility of topological codes implies that they 
can be defined in many lattices. This feature makes nat- 
ural the introduction of code deformations [H, [2(| [2^| , 
which are briefly described next. 

When two codes are very similar, in the sense that 
they differ only by a few qubits, it is possible to trans- 
form one onto another by manipulating these few qubits. 
This is specially natural for local codes that only differ 
locally. In particular, such local code deformations will 
not compromise the protection of the encoded informa- 
tion. These ideas where first explored in [l3| . where the 
geometry of a surface code is transformed step by step 
to compensate a change in the lattice geometry provoked 
by a transversal gate. In [l3[ code deformations are also 
used to initialize the code with an arbitrary state. This 
is done by 'growing' the code from a single qubit encod- 
ing the desired state. Notice that in this case encoded 
information is not protected on the early stages of the 
code deformation, when the code is still small. 

In general [20] , two main kinds of deformations may be 
distinguished: those which change the number of encoded 
qubits and those which do not. The formers can be used 
to initialize and measure encoded qubits, and the latter 
to perform operations on encoded qubits. In the case of 
topological codes, code deformations amount to changes 
in the geometry of a lattice, which may ultimately be 
understood as changes in the geometry of a manifold. 
When the topology of the manifold changes, initializa- 
tion or measurement of encoded qubits will happen in a 
well-defined way [2(| ■ When the manifold undergoes a 
continuous transformation that maps it to itself, a uni- 
tary operation is performed on the encoded qubits [20| . 
This unitary operation only depends on the isotopy class 
of the mapping. 

Code deformation can be naturally integrated with the 
successive rounds of stabilizer measurements mentioned 
in the previous section. In particular, as long as the 
deformations are local, one can perform them simply by 
changing the stabilizers to be measured at each stage of 
error detection [20j | . 

E. String operators 

In this and subsequent sections we only consider 2- 
dimensional topological codes with qubits as their basic 
elements, because the topological subsystem codes in- 
troduced in section IV fall into this category. For the 
same reason, subsystem code language, as opposed to 
subspace, will be used. 

In known 2D topological codes the logical Pauli opera- 
tors Xu . . . , Zk 6 N(£) can be chosen to be string opera- 
tors [HI, EH ■ These are operators O s with support along 
a set of qubits s that resembles a closed string. There 
are several types of strings, labeled as {k}. Two strings 
s, s' of the same type that enclose a given region, like a, b 
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FIG. 2: In this figure the geometry is that of a torus, with 
opposite sides of the square identified. The colored curves 
represent the support of string operators, with color stand- 
ing for strings labels. The strings a and b enclose a region 
and thus are homologically equivalent, producing equivalent 
operators. The string c encloses a region and thus is homolog- 
ically trivial, producing a stabilizer element. The strings d, e 
and / are homologically equivalent but have different labels, 
producing inequivalent operators. 

in Fig. [21 give equivalent operators O s O' s 6 S 1 . In other 
words, only the homology of the strings is relevant. In 
particular, boundary strings, those that enclose a region 
like c in Fig. [2j produce operators in S' . Moreover, S' is 
generated by boundary strings of some minimal regions. 
When two strings s, s' cross once, like a and d in Fig. [21 
O s and O s i commute or not depending only on the labels 
of the strings. Finally, two strings s, s' with a common 
homology class can be combined in a single string s" of 
a suitable type, in the sense that O s O S '0 S " € S' . For 
example, d and e in Fig. [5] can be combined in a string / 
of a suitable type. 

All this properties can be captured in a group L' ~ Vt, 
for some t that depens on the code. For example, t = 1 
for surface codes and t = 2 for color codes. The group 
L := L'/(il) has as its elements the string types {li} 
and the product corresponds to string combination. The 
comutation rules are recovered from L': crossing strings 
commute or not depending on whether their labels com- 
mute in L'. There are 2t labels l\, . . . , It that generate L. 
In a given manifold 2 — x nontrivial cycles that generate 
the homology group can be chosen. Each of them gives 
rise to 2t string operators with labels k, and the total 
2t(2 — x) string operators generate N(Q)/S' , so that the 
number of encoded qubits is k — t(2 — x)- 

Instead of strings, in general one can consider string- 
nets, where the strings meet at branching points [14j . 
The allowed branchings are those in which the product 
of all the labels involved is trivial. String-nets do not play 
a significant role in closed manifolds, but can be essential 
when the manifold has boundaries [14j . 

Let us check the criterium for topological codes of sec- 
tion MI Bl using the string operator structure. To this 
end a notion of connectedness is needed, but this can be 
obtained from the local generators of S': two regions or 
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sets of qubits are disconnected from each other if no local 
generator has support on both of them at the same time. 
Let Qo denote the support of an operator O. Then if 
O £ N(S) and Q = Qi U Q 2 with Q x and Q 2 dis- 
connected it follows that O — 0\0 2 and Q\ = Qon 
Q2 = Qo 2 f° r some 0\,0% G N(S). As for simply 
connectedness, it is easier to introduce a wider notion 
of 'trivial' region. A set of qubits Q forms a trivial re- 
gion when there exist string operators Oi that generate 
N{G)/S' and such that Q Qz n Q = 0. If O e N(S) is 
such that Qo is a trivial region then [O, Oi] — for the 
corresponding string generators Oi and thus O G Q. The 
criterium is satisfied within this language: if O € N{S) 
and Qo = QiU- ■ -UQ t with the Qi pairwise disconnected 
and each of the Qi a trivial region, then O G Q. 

F. Anyons 

Topological codes can be described in terms of string 
operators because they describe ground states of topo- 
logically ordered quantum models, that is, systems with 
emergent abelian anyons |l2| . Anyons are localized quasi- 
particles with unusual statistics. String operator repre- 
sent quasiparticle processes, and their commutation rules 
are directly related to the topological interactions of the 
anyons. Moreover, when an open-ended string operator 
is applied to the ground state, a pair of anyons is cre- 
ated on the ends of the string. The labels of the created 
anyons are those of the string, so that string labels are 
also quasiparticle labels. 

From the perspective of the code quasiparticles cor- 
respond to error syndromes, signaling a chain of errors 
along the string jl3j . Thus, keeping track of errors in 
a topological code, recall section UlI C[ amounts to keep 
track of the wordlines of these quasiparticles. Error cor- 
rection will success if the wordlines are correctly guessed 
up to homology [l3j . 

G. Boundaries 

From a practical perspective, codes that are local in 
closed 2-manifold like a torus are not very convenient. 
Instead, one would prefer to have planar codes. Thus, a 
way to create nontrivial topologies in the plane is needed, 
and this is exactly what is gained by introducing bound- 
aries. 

In a given code, different types of boundaries are 
possible. To start with, one can always consider ran- 
dom, structureless boundaries. The introduction of such 
boundaries will typically produce many local encoded 
qubits along the boundary. But these qubits are unpro- 
tected, and thus essentially useless. 

More interestingly, boundaries with well-defined prop- 
erties and non-local encoded qubits are also possible 
[l3l The defining property of such boundaries is 

that strings s with labels from certain subset M C L are 



^0 - 












Jb 










»J^f e 1 




c 1 





FIG. 3: This figure illustrates boundaries on 2D topological 
codes, which are displayed as dashed thick lines. The strings 
a and b, together with the boundaries, enclose a region and 
thus produce equivalent operators. The string c can end in the 
boundary because it can be decomposed in two strings that 
can end on it. The strings d and e enclose regions and thus 
produce stabilizer operators. The string / produces an stabi- 
lizer element or an undetectable error, depending on wether 
its label is allowed in the boundary. 



allowed to end in them, see Fig. [3j in the sense that O s 
belongs to the normalizer N{S). In other words, the in- 
troduction of the boundary changes the notion of closed 
string, by allowing on the boundary loose ends of strings 
of suitable types. The notion of boundary string also 
changes. Two strings s, s' of the same type that, to- 
gether with boundaries in which they can end, form the 
the boundary of a given region, as a and b in Fig. [3J 
produce equivalent operators so that O s O' s G S' . 

Notice that M should be a subgroup of L, because if 
strings with labels /, /' G M can end in the boundary then 
so can strings with label IV by splitting before reaching 
the bondary, as string c in Fig. [3J Also, any two la- 
bels I, V G M must commute in L'. In other case, the 
stabilizer would contain anticommuting elements, which 
is not possible. This is illustrated by the strings d and 
e of Fig. |31 which must produce commuting operators. 
Finally, L' should be maximal in the sense that for any 
I G" M there exist some V G M such that I and V an- 
ticommute in L' . In other case, according to the rules 
stated above, an Z-string s that surrounds a M-hole, like 
/ in Fig. [3j produces an operator O s that has to be- 
long to N{S) — S' because it is not a boundary but for 
which there is no other string s' such that {O s ,0' s } = 0. 
This is a contradiction. It is in fact possible to relax this 
last maximality condition, but at the cost of getting a 
boundary between two topological codes, rather than a 
boundary between a code and the 'vacuum'. 

Remarkably, boundaries in topological codes are di- 
rectly related to anyon condensation in the corresponding 
topologically ordered models [3Cj ■ It will become appar- 
ent in section llV El that this has important consequences, 
because only bosons can condense and this forbids cer- 
tain types of boundaries. 
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FIG. 4: (a) Part of a 2D lattice A with triangular faces and 
3-colorable vertices, (b) The lattice A derived from A. It is 
obtained by separating the triangles of A and adding one face 
per edge and vertex of A. Its edges are classified in three 
types. Solid edges are Z-edges, dashed edges are K-edges and 
dotted edges are X-edges. There is one qubit per vertex and 
the generators of the gauge group are related to edges. They 
are 2-local operators of the form XX, YY or ZZ, depending 
on the edge type. 



IV. A FAMILY OF TOPOLOGICAL 
SUBSYSTEM CODES 

The subsystem codes introduced in this section have 
their origin in a spin- 1/2 quantum model that shows 
topological order [3l| . The Hamiltonian of the model 
is a sum of 2-local Pauli terms that here will become the 
gauge generators. The Pauli symmetries of this Hamil- 
tonian where already described in 3l| , and thus to some 
extent the codes where already implicitly considered in 
that work. Here we explicitly work out all the details 
from the code perspective. In addition, diverse aspects 
that are important in practice are explored, such as the 
possibility of introducing boundaries and the computa- 
tional power of code deformations. 



A. Lattice and gauge group 

The family of codes Ca of interest is parametrized by 
tripartite triangulations A of closed 2-manifolds, not nec- 
essarily orientable. That is, each code Cv is obtained 
from a 2-dimensional lattice A such that (i) all faces 
/ £ F are triangular and (ii) the set of vertices V can 
be separated in three disjoint sets, in such a way that 
no edge e £ E connects two vertices from the same set. 
Fig. Ufa) shows an example. Alternatively, A is the dual 
of a 2-colex [32j . Following the notation used in previous 
works, the three sets of vertices are colored as red, green 
and blue. The faces of A will be simply called triangles. 

The first step in the construction of Ca is to derive 
a new lattice A from A, as exemplified in Fig. HJb) . In 
going from A to A, the triangles of A separate from each 
other giving rise to new faces. In particular, each of the 
edges and vertices of A contributes a face to A. The edges 
of A are divided into three subsets, E = Ex U Ey U Ez- 
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FIG. 5: (a) The four possible configurations at a given ver- 
tex for allowed subgraphs of A. A different Pauli operator 
corresponds to each of them, (b) A subgraph 7 (thick lines) 
of a lattice A (thin lines), obtained from a regular triangular 
lattice A (lightest lines). 



In Fig.|4jb), AT-edges are dotted, F-edges are dashed and 
Z-edges are solid. The Z-edges form the triangles of A. 
Each edge in A contributes a X-edge and a F-edge, in 
such a way that no two AT-edges or two F-edges meet. 
There are thus two ways to choose the sets of X and Y 
edges. 

The definition of Ca is now at hand. First, physical 
qubits correspond to the vertices of A. Second, the gauge 
group is Ga '■= (il) ■ (C e ) e6 £, with generators G e related 
to the edges e of A. These take the form G e := <J v cr v > 
for e £ E a , a = AT, Y, Z, where v, v' are the vertices con- 
nected by e. Thus, the generators are 2-local. This is an 
improvement with respect to previously known topolog- 
ical stabilizer codes, which have generators of weight at 
least 4. 



B. String operators 

This section describes N(Qa) and its center S' A = 
Qa H N{Qa) in terms of the string operator framework 
of section IIII El which is valid for these codes. It turns 
out in particular that V ~ V\, as in surface codes. In 
section HV El it will be apparent, however, that there ex- 
ist differences between the nature of the strings of Ca 
codes and those in surface codes. These are not captured 
by V ', which indeed does not contain all the information 
about the corresponding topological order. The details 
of the statements made in this section can be found in 
appendix [XJ 

We first seek a graphical representation of N(Ga)- 
Take any subgraph 7 of the graph of A, such as the one in 
Fig. Hlb), that has at each of its vertices one of the con- 
figurations of Fig. [H^a). This graph 7 produces a Pauli 
operator 7 = (g) v a v with a v = 1,X,Y,Z according 
to the correspondence of Fig. [5]Ja). Observe that such 
operators 7 belong to N(Ga)- Up to a phase, the cor- 
respondence between the elements of N(Qa) and graphs 
is one to one. 

These graphs either contain all the edges of a triangle 
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FIG. 6: A string operator as a subgraph 7 of A, displayed in 
thick lines. Its triangles come in pairs, each of them connect- 
ing vertices of the same color in A. 









'X^XXXX 


x^> 


SxSxSx 







FIG. 7: Two examples of string operators S£ related to ver- 
tices v in A. The string b shares color with the vertex that 
encloses, whereas for a the two colors differ producing a more 
involved operator. 



or none of them. Thus, each graph 7 determines a subset 
of triangles T 7 of the original lattice A. In Figs. [5][7] 
this subset appears shaded. Notice that the number of 
triangles of T 7 meeting at each vertex is even. In fact, 
any subset of triangles that meets this property can be 
realized as T 7 for some 7. 

String operators are obtained from string-like graphs 
such as the one in Fig. O Notice in the figure how tri- 
angles can be paired in a specific way. These pairs of 
triangles connect always vertices of the same color from 
the original lattice A. These allows to classify strings 
accordingly with labels r, g, b. It is a simple exercise 
to check that crossing strings operators commute if they 
have the same color and anticommute in other case, in 
accordance with L' . String-nets can be formed by intro- 
ducing branching points where three strings of different 
color meet. 

The group S' A is generated by small string operators 
related to vertices v of the original lattice A. In par- 



ticular, let us set <S° = 1 



r,g,b, with 7 the c- 



colored string going around v, as shown in Fig. [7J Then 
S' A = (il) ■ (S%) Vl c- These generators are only subject to 
the relations 



where the first product runs over the three colors and 
the second over the vertices of A. As a consequence, the 



rank of 5a is s = 2\V\ — 2. Since the number of encoded 
qubits is k = 2 — x, it follows that the number of gauge 
qubits is r = n - k - s = \3F\ - 2\V\ + X = 2|F| - x, 
showing that gauge qubits see the global structure of the 
manifold. 

What can string operators say about the code distance 
dl Given an operator 1 £ N(Q), consider the subset E' 
of the edges of 7 with elements all the X- and F-edges of 
7 and one of the three Z-edges that correspond to each 
triangle in T 7 . Then G — YieeL £ 6 a and it is easy 
to check that |0 7 G| = |T 7 |. Therefore d < dx, with cZt 
the minimal length, in terms of the number of triangles, 
among nontrivial closed strings. A lower bound for d is 
given in the next section. 



C. Homology of errors 

This section offers a homological description of error 
correction for C\. The main idea is that the error syn- 
drome can be identified with the boundary of errors, con- 
sidered as paths on the surface. Then error correction 
succeeds if this paths can be guessed up to homology. 
It is worth noting that the notation and results in this 
section will not be used again. 

To fix notation, we recall first some basic notions. Let 
A denote the additive group of Z2 1-chains in A. Its 
elements are sets of edges 6 C E and addition is given by 

5 + 6' = (6 U 5') - (6 n 5'). The boundary 06 of S e A 
is the set of vertices in which an odd number of edges 
from 6 meet. The elements 5 £ A with 85 = are called 
cycles and form a subgroup Z c A. Boundaries form 
a subgroup B C Z, generated by elements of the form 

6 = {ex, e2,ez} with the three edges of a given triangle. 
The first Z 2 homology group of A is H x := Z/B ~ Z§ 
with h = 2 — x the number of independent nontrivial 
cycles of the closed surface formed by A. Two chains 
<5, 6' £ A are said to be equivalent up to homology, 5 ~ 6', 
if 5 + 5' £ B. 

Consider a morphism / r : V n — > A defined by f r (i) — 
and the following action on single qubit operators A c , 
Yv, where V £ V. A anticommutes exactly with two 
operators of the form S*, v £ V. The corresponding two 
vertices are connected by an edge e £ E, and / r (A c ) = 
{e}. f T (Yv) is defined analogously. It is easy to check 
that /r[£?A] = B and that for any O £ "P n the set df T (0) 
contains those vertices v £ V such that {O, S^} = 0. 
Moreover, if 7 is a string then / r (0 7 ) £ Z and if 7 is red 
/ r (0 7 ) £ B. Indeed, if {7^} is the set of red strings, then 
/ r gives an isomorphism N(Sa)/(Ga • (0 7i ),-) ~ Hi, 

Consider in addition an analogous morphism /b with 
blue color playing the same role as red in / r . Then for 
any O £ V n we have O £ 6a if and only if f c (0) £ B for 
c = r, b. Similarly, O £ N(S A ) if and only if f c (0) £ Z 
for c = r, b. This shows that error correction will succeed 
as long as errors can be guessed up to homology. In 
detail, suppose that the code suffers a Pauli error O. The 
error syndrome can be expressed in terms of the two sets 
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FIG. 8: The proposed ordering for the measurements of the 
edge operators. It does not depend on the particular geometry 
of the lattice A because it is dictated by the coloring of its 
vertices. 



df T (0),dfh(0) C V. Suppose that an attempt is made 
to correct the errors by applying some O' G V n such that 
dfc(O) — df c (0') for c = r,b. Then error correction 
succeeds if and only if O'O G Q, that is, if and only if 
/ c (0)~/ c (0') forc = r,b. 

Although error correction can be expressed in this ho- 
mological terms, this is really not the most natural thing 
to do, basically because it involves an arbitrary choice of 
two of the three available colors. In this regard, notice 
that not any set of edges S E A can be obtained from an 
operator O £ V n as S — f T (0), and that the cardinalities 
of / r (0) and fb(0) by no means are enough to compute 
\0\. This makes unfeasible a direct translation of the 
ideas used in |l3[ for error correction in surface codes. 

In order to give a lower bound for the distance d of 
the code, the definition of the mappings / c , c = r,b 
must be modified. We set /c(0 c cr c ) := J2 V fc( a v), 
where <7 C = 1 B , AT C , Yg, Zy, and we fix f' c (lv) '■= 0, 
f&Xv) := f c (X v ), f c {Y v ) := / c (Y fl ) and f c (Z,) is de- 
fined in analogy with f c (Xy). The new mappings f c are 
not group morphisms, but they do keep the good prop- 
erties of the f c mappings listed above. And they satisfy 
\0\ > |/c(0)|, which immediately leads to the bound 
d> <1l with dh the minimal length, in terms of the num- 
ber of edges, among nontrivial closed loops in A. 



D. Syndrome extraction 

As indicated in section fill C[ in a topological quantum 
memory one has to keep track of errors by performing 
round after round of syndrome extraction. This raises the 
question of how fast and simply the stabilizer generators 
of a code Ca can be measured. The faster the measure- 
ments the less errors will be produced in the meanwhile, 
and the simpler they are the less faulty gates they will 
involve. Of course, what fast and simple really mean will 
depend on particular implementations, that is, in the ba- 
sic operations at our disposal. 

To keep the discussion general, take gauge generator 



FIG. 9: A hypothetical geometry for a Ca code with three 
holes of different colors. According to the properties of bound- 
aries, the string operators a and b are equivalent up to sta- 
bilizer elements. This is a contradiction because due to the 
way they cross they anticommute. 

measurements to be the basic components of the syn- 
drome extraction. At each time step the measurement of 
any subset of generators {Gi} is allowed as long as each 
physical qubit only appears in one of the C;. Then, in 
any code C\ it is possible to cyclically measure all the 
stabilizer generators by performing six rounds of mea- 
surements. The time step at which each generator is to 
be measured is indicated in Fig. [5J Notice that Z-edges 
are measured at even times and X- and y-edges at odd 
times. From the time steps 1-3 the eigenvalue for opera- 
tors Sy at blue vertices are obtained, from the steps 3-5 
those for red vertices, and from the steps 5, 6 and 1 (this 
last one in the subsequent cycle) those for green vertices. 
It is not clear whether this number of time steps is op- 
timal, since in principle 4 or 5 could be enough. As a 
comparison, the number of steps needed for Bacon-Shor 
codes is 4. In this sense the 6 steps are not bad, tak- 
ing into account that the codes Cy do not benefit from 
the the separation of gauge and stabilizer generators into 
X-type and Z-type, as Bacon-Shor codes do. 



E. The problem of boundaries 

This section shows why it is not possible to introduce 
boundaries with the properties discussed in section UlI GI 
This has important practical consequences, since there is 
no other known way to introduce a nontrivial topology 
in a completely planar code. Notice however that we can 
always flatten a manifold to get a 'planar' code, at the 
price of doubling the density of physical qubits in the sur- 
face. Also, the absence of boundaries makes less practical 
the use of code deformations, although they are still pos- 
sible, as shown in section HV Fl In any case, this leaves 
open the question of whether other kind of interesting 
boundaries can be introduced. 

The existence of boundaries leads to the following con- 
tradiction. According to the properties listed in section 
IIII GI there are three potential kinds of boundaries, one 
per color. Each of them only allows strings of its color 
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to end on it. Clearly either all the boundaries can be 
constructed or none of them. Thus, suppose that the 
three of them are allowed and consider a geometry as the 
one in Fig. [SJ with three holes, one of each color. Take 
a string-net 7 that connects the three holes, as in the 
figure, and deform it to another string-net 7'. It follows 
from the properties of string operators and boundaries 
that 7 7 ' € S', but also that {0 7 , Oy} = since they 
cross at a single point where they have different color. 
This is not possible. 

In section ITVB I it was noted that the string label group 
V is the same in surface codes and the subsystem codes. 
Since according to section Ull Gl the set of allowed bound- 
aries is dictated by L', it could be expected that surface 
codes would neither have boundaries. However, this is 
not the case: two kind of boundaries can be constructed 
in surface codes The point is that there is a key 

difference between the two families of codes: in surface 
codes the three types of strings are not equivalent in any 
sense, so that the previous reasoning is not valid. 

At a deeper level, this difference between the codes 
has its origin in the difference between the correspond- 
ing topological orders. Indeed, L' does not encode all the 
information about the properties of anyons. In surface 
codes two of the quasiparticle types are bosons, and the 
third a fermion [12| . whereas in the subsystem codes the 
three are fermions [31j . The connection between anyon 
condensation and boundaries is thus crucial: nice bound- 
aries cannot be introduced in these topological subsys- 
tem codes because all the string operators are related to 
fermions, which cannot condense. 



F. Code deformation 

This section explores the potential of code deforma- 
tions in the topological subsystem codes C\. We show 
how initializations and measurements of individual logi- 
cal qubits in the X and Z basis are possible through cer- 
tain topology-changing processes on the manifold. And 
also, that CNot and Hadamard gates can be in princi- 
ple implemented through continuous deformations of the 
manifold, but not in a practical way. 

To begin with, a manifold and a set of logical operators 
Xi, Z\, ■ ■ ■ X/., Zf. must be selected. We choose a /i-torus, 
that is, a sphere with h holes. Codes C\ on such a man- 
ifold provide 2h logical qubits, but only h of them will 
be used, with the choice of logical operators indicated 
in Fig. [Tola). The rest of logical qubits are considered 
gauge qubits. 

When a new handle is introduced, a logical qubit is 
created and initialized in a definite way, see Fig. [TUJb). 
In the figure the two surfaces are supposed to be already 
connected, so that a handle is really created. There are 
two ways to introduce a new handle in a surface such as 
that of Fig. [TUTa). depending on whether the process of 
Fig. ITUTb) occurs inside' or 'outside' the surface. In the 
former case the new qubit is initialized in a Z eigenstate 




FIG. 10: (a) A sphere with h holes can encode 2h qubits, 
but we choose to encode just h. The logical Xi, Zi operators 
correspond to the strings in the figure. Red strings give X's 
and blue strings Z's. (b) When the topology of the surface 
changes as indicated here two qubits are introduced in the 
code. They are initialized in a fixed way. In particular, the 
string operator in the figure is a boundary before the deforma- 
tion takes place and thus has a fixed value. This value is not 
changed by the deformation because it occurs in a different 
part of the code. 



and in the latter in a X eigenstate. Whether the initial- 
ization occurs in the Z or X basis depends on which of 
the two string operators of the new qubit was a boundary 
initially. This operator has its eigenvalue fixed before the 
deformation occurs, and during the process it is topolog- 
ically protected at all times 20]. The particular sign of 
the eigenstate depends on the arbitrary sign choices for 
the logical operators and S\. If the initialization process 
is reverted it yields a measurement in the corresponding 
basis [H. 

It is always possible to detach a qubit, a torus, from 
the rest of the code. This does not involves any mea- 
surement, because the strings running along the cutting 
line are boundaries [20]. Similarly, there is no problem 
in attaching a torus to the code to add a logical qubit. 
But once a logical qubit is isolated, it can undergo code 
deformations independently. Consider a mapping that 
exchanges the two principal cycles of the torus and shifts 
the lattice a bit if necessary to adjust the color correspon- 
dence, for example by rotating the torus. Such a map- 
ping can exchange X and Z operators, which amounts to 
a Hadamard gate. There exist an important drawback, 
though. This deformation cannot be realized in 3D with- 
out producing self-intersections of the surface. Still, it 
is conceptually interesting that the Hadamard gate can 
be obtained from purely geometric code deformations be- 
cause this is not possible in surface or color codes, where 
X and Z type operators correspond to different type of 
string operators and a transversal Hadamard gate must 
be added to the picture [l3j. Because color is just a mat- 
ter of location in the lattice, strings of different colors are 
equivalent up to lattice translations. This is in essence 
what makes it possible the geometric implementation of 
the Hadamard gate. 

A controlled phase gate 1 — |(1 — Zi)(l — Zj) on a 
pair of logical qubits i,j can be implemented through 
a 'continuous' deformation of the code. The process is 
indicated in Fig. 1111 It follows from the way in which 
the logical operatorators evolve [2(| that the complete 
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FIG. 11: The deformation that produces a controlled phase 
gate. (a)The code before the deformation takes place and 
a particular string, (b) The deformation moves one of the 
'holes' in the top part around the other, as indicate by the 
solid line with an arrow. To recover the original shape, as 
indicated by the dashed line, the two 'tubes' have to overlap 
unavoidably, (c) After the deformation, the string operator 
has been mapped to the product of these two string operators. 



process amount to a controlled phase gate, up to some 
signs in the final logical operators that depend on the 
choice of 5a. A CNot gate can then be obtained by 
composing this gate with Hadamard gates. But again, 
such a code deformation requires to overlap the surface 
of the code with itself, see Fig. [TTT b) . 



V. STATISTICAL PHYSICS OF ERROR 
CORRECTION 

In [l3T | an interesting connection between error correc- 
tion thresholds for surface codes and phase transitions in 
2D random bond Ising models was developed. Similar 
mappings exist also for color codes 22], in this case to 
2D random 3-body Ising models. In both cases, the CSS 
structure of these codes is an important ingredient in the 
constructions: they are subspace codes with S = SxSz 
in such a way that S a is generated by products of a op- 
erators, a = X, Z. To take full profit of this, the noise 
channel for each qubit must be a composition of a bit- 
flip channel £bf(j>) '■= {(1 — p) 1 !,? 1 X} and a phase flip 
channel £ p {(p) ■= {(1 — p)^l,p^ Z}. 

There are two main obstacles to construct a similar 
mapping for the codes Ca- The first is that they are sub- 
system codes rather than subspace codes. The second, 
that the gauge group cannot be separated in a A and a 
Z part. As we show below, both can be overcome. 



A. Mapping to a statistical model 

Rather than directly considering the codes Ca, this 
section deals with the general mapping from any given 
stabilizer subsystem code to a suitable classical statisti- 
cal model. For simplicity, each qubit in the code is sup- 
posed to be subject to a depolarizing channel Sdcpip) '■= 
{(1 - p)H, (p/3)ix, (p/3)*Y, (p/3)sZ}, with p the er- 
ror probability, but more general channels are possible 
within the same framework. 

To build the classical Hamiltonian model, the first step 



is the choice of a set of generators {Gi}\ =1 of Q/(il). 
These generators can be captured in a collection of num- 
bers gfj = 0, 1 defined by 



Gi<7j = : !r">,C, 



(6) 



with a = X,Y, Z, i — 1, . . . , I, j = 1, . . . , n. Attach a 
classical Ising spin = ±1 to each of the generators Gi . 
The family of Hamiltonians of interest is 



H r (s):=-J E E^/II 

a=X,Y,Z j=l i=l 



(7) 



with parameters tJ — ±1 such that T f T j T j — 1- The 
coupling J > is introduced to follow conventions. No- 
tice that codes with local gauge generators give rise to 
local Hamiltonian models. Since gfj+gYj+gfj = m °d 2 
the Hamiltonian can be rewritten as 

H r (s) = n EC 1 + T ? II ^ X 1 + T J II "ih ( 8 ) 



The partition function for these Hamiltonians is 

Z(K,t) = ^~V^ (s) (9) 

with K := (3 J and j3 the inverse temperature. 

The goal is to express the class probabilities p(E), E E 
V n , in terms of the partition function ([9]) for a suitable 
t. Let t = te be such that 



£<x 0i— y—. (io) 

3 

Similarly, for each G G Q choose any s — sq such that 



GaGj 2 ■G 



(11) 



We write s" = s's if s" = s'^Sj. Then sg s G' = sgc and 
for any spin configuration s, E € V n and G € Q it can 
be checked that 



h t eg {s) = H TE (s G s) 



(12) 



In the depolarizing channel the probability for a Pauli 
error E is p(E) = (p/3)' s '(l —p) n ~\ E \. It may be written 
as 

p(E) = c -« e -^ B (^) ; (13) 
where c p := e 3Kp + 3e~ Kp and /3 P :— K p /J with 

1 ~P 

The desired connection follows from (fT2|) and (p~3|) , which 
give 

p(E) - E P^ EG ) = ^ Z ( K ?> T£ )' ( 15 ) 



Gee 



where w is the number of redundant generators of Q, that 
is, w = I — I' with I' the rank of Q / (il). 
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B. CSS-like codes 

To connect the results of the previous section with the 
works [HI, H3, CSS codes must be considered. These 
are codes with Q = {il)QxQz f° r some Q a generated by 
products of a operators, a — X, Z. And, instead of a 
depolarizing channel, the noisy channel must take the 
form 8 — £bf{p) ° £pf(p')- This allows to treat X and Z 
errors independently [l3|. Here we consider the case of 
bit flip errors, phase flip errors are analogous. 

The construction is similar to the one in the previous 
section. It starts with the choice of generators Qx = 
(Gi) l i=1 . The relevant Hamiltonians read 



n I 
II': -^-lK- 

3=1 i=l 



(16) 



where Tj := and := gfj. The probability of an 
error E that is a product of X operators is 



— Z(K' p ,T E ), (17) 



where w is the number of redundant generators of Qx 
and 



-2K' . V 



1-p 



(18) 



defines the Nishimori temperature [26j. 

Bacon-Shor codes provides an example of gauge CSS- 
like codes. With the above procedure, they yield models 
that amounts to several copies of the ID Ising model. 



C. Symmetries 

Interestingly, the redundancy of the generators of Q is 
directly connected to the symmetries of the Hamiltoni- 
ans ([7]). Suppose that the generators are subject to a 
constraint of the form 



for some set of indices I. Then (TT2"]) gives 
H T (s) = H T (sl[s Gi ). 



(19) 



(20) 



iei 



In other words, making the most natural choice for the 
sg i it follows that the Hamiltonian is invariant under the 
transformation 



Si 5 



i e I, 
i (£1. 



(21) 



Thus, global constraints lead to global symmetries and 
local constraints to local symmetries. 



As a particular example, consider surface codes, which 
are mapped to Ising models [13j . In these codes the prod- 
uct of all X-type stabilizers equals the identity, producing 
a symmetry that is simply the global Z2 symmetry of the 
Ising model. 



D. Error correction and free energy 

We now put equation (|15[) to use in the error correc- 
tion framework of section QTEj Recall that, after the syn- 
drome has been measured, one has to find the most prob- 
able class of errors among several candidates Ei := EDi. 
This amounts to compare the probabilities p(Ei) or, al- 
ternatively, the quantities Z(K p , trj). And to do this, it 
is enough to know the free energy differences [l3| 

\(K p ,r E ) := pF(K p ,r Ei ) - (3F(K p ,t e ), (22) 

where F(K,t) — — T log Z(K, r) is the free energy of a 
given interaction configuration r. For example, in the 
Ising models that appear for 2D surface and color codes 
these are domain wall free energies. 

In practice, the computation of (I22|_mav be difficult. 
In this regard, it has been suggested [13l |. in the context 
of surface codes, that in the absence of glassy behavior 
the computation of (|2"2"|) should be manageable, and in 
[2lT | a possible approach was sketched. 



E. Error threshold and phase transition 

In surfaces codes there exists an error probability p c , 
the error threshold, such that the asymptotic value of 
the success probability po, in the limit of large code in- 
stances, is one for p < p c and l/4 fc for p > p c [21]. This 
is directly connected to an order-disorder phase transi- 
tion in a model with random interactions. An analogous 
transition is observed for the random model that corre- 
sponds to color codes [13, [H, HH . It is then natural to 
expect a similar connection in other topological codes, as 
we describe next. 

Consider a random statistical model with Hamiltonian 
([7]) in which the parameter r is a quenched random vari- 
able. That is, r is random but not subject to thermal 
fluctuations. The probability distribution p(r) is such 
that the signs of rf and tJ are independent if i ^ j. 
For each i, the case rf- = = 1 has probability 1 — p 
and the other cases have probability p/3 each. In other 
words, if r = t e then p[r) = p{E) with p(E) given by 
the depolarizing channel £dcp(p)- 

In thermal equilibrium the model has two parameters, 
the temperature T and the probability p. For the map- 
ping only a particular line in the p-T plane is relevant, 
the Nishimori line [26], given by the condition K = K p 
that has its origin in (TT3")) . The error correction success 
probability in can be written in terms of this statis- 
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tical model as follows: 



Pa 



4k 



-Ai(K p ,r) 



i=2 



(23) 



K„ 



where [-]k p ■— ^2 T p( T )' denotes the average over the 
quenched variables. 

Suppose that the code has a threshold probabilit y p c 
below which p' 1 in the limit of large codes. Then [35| , 
in the random model the average of the free energy differ- 
ence (f22j) diverges with the system size, [Ai(K,T)]x„ —> 
oo, f or p < p c along the Nishimori line. This is exempli- 
fied [1 31 ] by surface codes and the corresponding random 
2d Ising models, where [Ai(K, t)]k p is the domain wall 
free energy. It diverges with the system size below p = p c 
and attains some finite limit over the threshold, signaling 
an order-disorder phase transition at p c . A similar be- 
haviour can be expected for other topological codes. For 
2D color codes this was shown in 1231 ■ 



F. The Hamiltonian model for Ca codes 

The above mapping can be immediately applied to the 
subsystem codes Ca- Choose as generators of the gauge 
group the edge operators O e , so that there is an Ising 
spin s e at each edge e. The Hamiltonian takes the form 

n 

H^(s) := -J^lf 828384. + T? 818384, + T?8!8 2 , (24) 

i=i 

where the sum runs over vertices and for each of them 
the Ising spins s±, S2, S3, S4 correspond respectively to the 
X, Y and two Z edges meeting at the vertex. 

The Hamiltonian (IM1) has a local symmetry at each 
triangle. In particular, flipping the three Ising spins of 
the triangle leaves invariant. This is so becuase the 
product of the three edge operators in the triangle equals 
the identity. There exist also a global Z2 x Z2 symmetry 
that follows from the global constraints in ([5]). The local 
constraints in ([5]) do not provide any symmetry as they 
are trivial in terms of the gauge generators. 



G. Faulty measurements 

The mapping considered up to know is only suitable 
if perfectly accurate quantum computations are allowed 
in error correction. This section generalizes it to include 
errors in the measurements of the stabilizer generators. 

Following 13], take as a goal the 'indefinite' preser- 
vation of the content of a quantum memory. Time is 
divided in discrete steps. At each time step, the memory 
suffers errors and at the same time the stabilizer genera- 
tors are imperfectly measured. Then if from the history 
of measurements one can correctly infer the actual his- 
tory of errors, up to a suitable equivalence, the memory 
is safe. 



The results in [13j, l2l| show that for surface codes there 
exists of a noise threshold below which long time storage 
is possible for sufficiently large codes. The same behavior 
can be expected for other topological codes, but the con- 
struction of a suitable random statistical model for each 
code is required first. Here we generalize the construction 
of [l3j j to subsystem codes and depolarizing channels. 



1. Depolarizing channel 

Consider first the case of a depolarizing channel £dcp(p) 
occurring for each physical qubit between each round of 
measurements. We adopt the convention that at a given 
time t first errors occur and then faulty measurements 
are performed. 

Recall that in the mapping of error correction to a sta- 
tistical model errors were mapped to interactions through 
the rj 7 , see (fit))) . The new elements here are time and 
faulty measurements. Since errors can occur at different 
time steps t, a time label must be added to the rj's to 



get the collection of signs r = (tZ), subject as before to 
the constraints t^t^t^. = 1. To represent errors in the 
measurements of stabilizers, first a set {Sk}™ =1 of gener- 
ators of S to be measured at each time step t must be 
chosen. Attach to them a collection of signs Kkt = ±1. 
The correct (wrong) measurement of the i-generator at 
time t corresponds to Kkt — 1 ( K kt — — 1)- In the statis- 
tical model the r and k are quenched variables, r follows 
the same distribution as before, dependent on the prob- 
ability p, and each Kkt is independent an takes value —1 
with probability q. For this to make sense under the 
mapping, errors in the measurements must occur inde- 
pendently with a fixed probability q. This will not be 
true in most settings. Still, it is a useful assumption be- 
cause knowing the correlations between errors can only 
improve error correction. In analogy with the gfj defined 
above, the stabilizer generators are captured in a collec- 
tion of numbers =0,1 defined by 



(25) 



with a = X, Y, Z, k = 1, . . . , m, j = 1, . . . , n. 

Recall also that in the original mapping gauge genera- 
tors Gi were mapped to Ising spins Sj . The reason for this 
was that gauge generators play the role of basic equiva- 
lences between errors. Now for errors that occur at the 
same time t these kind of equivalence happens again, rep- 
resented by spins sn- But in addition there is an equiva- 
lence between errors that involves errors at different time 
and measurement errors. If at times t an t + 1 a given 
error occurs and the measurements at time t of the sta- 
bilizers that would detect the error fail, then these errors 
altogether go unnoticed but produce no harm. Thus, 
two collections of errors that differ only by such an event 
should be considered equivalent. Therefore, Ising spins 
that represent this equivalence are necessary. This can 
be achieved by attaching two Ising spins s^ t , sj t to the 
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t-th time step and j-th qubit. The Hamiltonians are 
H T , K (s) := - ^EEE ^ < II" • 

a j t i 

-*ee*« n n . (26) 

where s^ t := s^sjt and the range of values of the different 
indices should be clear from the context. In order to 
recover the probability of a given set of errors from the 
partition function the relations 



3e -W = _£_ e -im _ _J_ 
1 — p 



(27) 



must hold. 

For each time step t, the Hamiltonians (|26|) keep the 
symmetries (|2"TT) . In addition, there is a symmetry for 
each gauge generator and time i'. Namely, 



5 i* 



Sit 



= j(-l)^- s j i , f=f, 

-s it , i = i', t = t',t' + 1, 



otherwise. 



(28) 



Therefore, local gauge generators give rise to a (random) 
gauge model. 



2. Bit flip channel 



VI. CONCLUSIONS 

Topological codes are intrinsically local, and gauge or 
subsystem codes can have interesting locality properties. 
In this paper we have introduced a family of topological 
subsystem codes, thus putting together the two concepts. 
The gauge group of the code is generated by 2-local op- 
erators, which compares well with surface or color codes 
that have at least 4-local and 6-local generators, respec- 
tively. In particular, the measurement of these 2-local 
operators is enough to recover the error syndrome. 

We have argued that these codes do not allow the in- 
troduction of boundaries with nice properties, which mo- 
tivates further research. There are probably interesting 
topological codes still to be discovered. One could look 
for example for subsystem codes with nice boundaries or 
with interesting transversality properties as those found 
in color codes. 

We have also explored a general connection between er- 
ror correction in subsystem codes and statistical physics. 
The connection is specially meaningful in the case of 
topological codes, where the error threshold maps to a 
phase transition in the corresponding statistical model. 
There is a lot of work to do in this direction. For exam- 
ple, the computation, probably numerically, of the error 
threshold of the topological subsystem codes presented 
here. 



Finally, consider the simpler case of a bit flip channel 
£bf{p) in a CSS-like code. As noted above, the case of 
a phase flip channel is analogous and if both channels 
happen consecutively they can be treated independently. 

The construction is an extension of the one in section 
IV Bl The Tj's and Si are respectively replaced by the 
signs Tj t and the Ising spins su- Given a choice of gener- 
ators Sz = (Sk) to be measured at each time step, there 
is a corresponding collection of signs rejy. The genera- 
tors take the form Sk = ± ®j Z^ h] for some hkj =0,1. 
There is also an Ising spin Sj t for each physical qubit j 
and time step t. The Hamiltonians read 



H'Js) := - JEE T * °j* II s u + 

j t i 



k t 



Instead of (|27|) . the right conditions are now 

-2/3K _ g 

~l-q 



-2/3J _ P 

1-p 

The analog of (J55J) is 



(30) 



s jt 



Sit 



l(-l)^ SJt , t = t>, 

[Sjt, t^t', 

-su, i = i',t = t',t' + 1, 
su, otherwise. 



(31) 
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Appendix A: Structure of N(Q) 



This appendix is a complement to section IIVB1 and 
uses the same notation. 

A color code [14[ can be obtained from a lattice 
A with the properties enumerated in section IIV Al The 
construction is the following. First, there is one qubit per 
triangle, so that the relevant Pauli group is V\f\- Given 
a collection of triangles T ~ set Xt '■— &>iX Ti , 

Zf := ® i Z Ti . If each vertex v £ V is identified with the 
set of triangles meeting at v, the stabilizer for the color 
code is := SxSz with Sx '■= (X v ) v , Sz ■— (Z v ) v . Let 
{Ti} be the collection of those sets of triangles that have 
an even number of triangles meeting at each vertex. Then 
N(S%) = {H)N X N Z with N x := (X Tt ) t , N z := (Zr t ) t -. 



16 



Next, consider the morphism / : N(Q A ) — > Nx such 
that f(0) = X Ty for any subgraph 7 of A and O £ N(Q A ) 
such that O a 7 . The kernel of / is formed by those 
operators that only involve Z-s, not X-a or Y-s. That 
is, ker/ = (il}(0^, v ) v C S\ where c v is the color of v in 
A. Since f{0^) = X v for c ^ c v , f[S' A ] — Sx and thus 
S' A I ker / ~ Sx ■ This implies that there are no other 
constraints for the generators of S' A apart from the ones 
in ([5]), because exactly two of the generators {X v } of Sx 



are unnecessary [ll]. Finally, it is easy to check that 
for any string operator X T £ N x , as described in flij . 
there exists a string-like graph 7 such that /(0 7 ) = Xt, 
so that / is onto and N(Q A )/S' A ~ Nx/Sx- Then the 
properties of the string operators in Ca are consequences 
of those for string operators in C A . This is in particular 
true regarding the generating set and the composition 
rules, but not for commutation rules, which have to be 
worked out separately. 
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